Phase diagram of a frustrated mixed-spin ladder with diagonal exchange bonds 
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Using exact numerical diagonalization and the conformal field theory approach,, we study the 
effect of magnetic frustrations due to diagonal exchange bonds in a system of two coupled mixed- 
spin (1, i) Heisenberg chains. It is established that relatively moderate frustrations are able to 
destroy the ferrimagnetic state and to stabilize the critical spin-liquid phase typical for half-integer- 
spin antiferromagnetic Heisenberg chains. Both phases are separated by a narrow but finite region 
occupied by a critical partially-polarized ferromagnetic phase. 



PACS numbers: 75.10.Jm, 75.10.Pq 

I. INTRODUCTION 

Many experiments on bimetallic quasi-one-dimensional 
(ID) molecular magnets imply that the magnetic prop- 
erties of these compounds are basically described by the 
Heisenberg spin model with antiferromagnetic exchange 
couplings .iS In the last decade, there has been an in- 
creasing experimental and theoretical interest in these 
mixed-spin systems exhibiting intriguing quantum spin 
phases and thermodynamic properties. In particular, a 
number of recent studies has been focused on ground- 
state properties of mixed-spin ladders^ 

In this paper we study the ground-state phase diagram 
of a ferrimagnetic two-leg ladder containing frustrating 
antiferromagnetic diagonal exchange bonds (see Fig. QJ. 
The model is defined by the Hamiltonian 

L 

^ = ^ [Jl (si„ • s 2n+1 + s 2 „ • Sl n+1 ) + J±s ln ■ s 2 „] 

n=l 

L 

+ J 2 (Si„ • Si„ +1 + s 2 „ • S 2 „+i) , (1) 
n=l 

where L is the number of rungs and the spin operators Si n 
and s 2n are defined on the rung with index n: (si„) 2 = 
si(si + l), (s 2n ) = s 2 (s 2 +l), si > s 2 , and J 1; J 2 , J± > 0. 
We introduce the frustration parameter a = J2/J1, the 
spin ratio a = si/s 2 , and set the energy and length scales 
by Ji = 1 and do = 1, where ao is the spacing between 
neighboring rungs. In the remainder of this paper, if not 
especially noted, we will consider the case Jj_ = J\. 

As a function of the frustration parameter a, the 
classical phase diagram of exhibits three phases de- 
scribed by the angles (9, 4>) fixing the directions of the 
classical spins Sx n and s 2n in respect to the classical 
ferrimagnetic configuration with spins si n and s 2 „ ori- 
ented along the axes Zi and z 2 , respectively (Fig. 
The classical canted state (C) shown in Fig. ^ is sta- 
ble in the interval a c \ < a < a c2 where a c \ = 3[(tr 2 + 
1) - y/{a 2 + l) 2 -32ct 2 /9]/8ct and a c2 = [-(a 2 + 1) + 
yj (a 2 + l) 2 + 32er 2 ]/8cr are second-order phase transi- 
tion points separating the C phase from the ferrimag- 



netic (F) (9, <j)) = (0, 0) and the mixed-spin collincar 
(9, 4>) = (7r/2,7r/2) phases, respectively. In the case of 
special interest (si,s 2 ) — (1, \), the classical transition 
points are a c \ = 0.3219 and a c2 = 0.4606. Notice that 
the square-lattice J\ — J 2 mixed-spin Heisenberg model 4 
exhibits similar classical magnetic phases which persist 
in the quantum phase diagram. On the other hand, the 
following analysis of Q implies that in one space di- 
mension the ferrimagnetic phase continues to exist in the 
quantum phase diagram, whereas the classical collincar 
magnetic state is completely destroyed by quantum fluc- 
tuations. Instead, for a > a 2c there appears a singlet 
quantum paramagnetic phase which is critical (gapped) 
for half-integer (integer) rung spins S« = Si„+s 2 „. As to 
the classical canted phase, it is argued that the longitu- 
dinal ferromagnetic order survives quantum fluctuations. 
On the other hand, on general grounds we may expect 
that the transverse magnetic ordering does not survive 
quantum fluctuations. A closely related phase diagram 
appears in a special class of lattice models with quantum- 
rotor degrees of freedom^ 
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FIG. 1: Sketch of the classical canted state for J± = Ji 
described by the angles (±#, i(f>) for every magnetic cell com- 
posed of two neighboring rungs. 8 and <j> measure the devia- 
tions of the classical spins from the ferrimagnetic configura- 
tion: si n ||zi and S2 n [|z2. 



The following analysis of the quantum phase diagram 
is performed by using exact numerical diagonalization 
(ED) of small periodic systems, finite-size analysis of the 
ED data, and analytical spin- wave calculations. The em- 
phasis is on the properties of the quantum paramagnetic 
phase. 
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II. MAGNETIC PHASES 

As may be expected, the positions of the classical phase 
transition points a c \ and a c2 are changed by quantum 
spin fluctuations. Using ED and a simple finite-size scal- 
ing, it is possible to find precise estimates for the quan- 
tum transition points. The latter are connected with 
the following changes in the total spin of the ground 
state St- St = (si — s 2 )L for < a < a c \ (F phase), 
< St < (si — s 2 )L for a cl < a < a c2 (C phase), 
and St — for a > a c2 (quantum paramagnetic phase). 
The extrapolated data for L = 8, 10, 12, and 14 give the 
results a c i = 0.341 and a c2 = 0.399 showing that the re- 
gion occupied by the quantum C phase is narrowed but 
definitely finite. Figure [3 provides a summary of the re- 
ported results in terms of the net ferromagnetic moment 
per rung Mq for L — 12. 
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FIG. 2: ED results for the ferromagnetic moment per rung 
Mo vs J2/J1 (dashed line). The step-like form of Mo(J2/Ji) 
is connected to finite-size effects. The midpoints of the steps 
close to a C 2 are well approximated by the ansatz Mo = 
m\{a C 2 — a) 13 + m2(« C 2 — a), where mi = 1.65, 1712 = 0.56, 
and P = I (solid line). 

The ferrimagnetic phase has already been studied 
for the model JQ) without frustration^ As Mq > 0, 
both magnetic phases (F and C) are characterized by 
quadratic spin-wave excitations 



E(k) 



M 



k 2 + 0(k 4 ) , 



(2) 



where p s is the ferromagnetic spin-stiffness constant i In 
approaching a c \ from the ferrimagnetic phase, the lin- 
ear spin-wave theory predicts that the lower spin wave 
branch softens in the vicinity of k = 7r and the gap 
at this point vanishes for a > a c \. Thus the linear 
Goldstone mode, characteristic of the classical C phase, 
seems to survive quantum fluctuations although on gen- 
eral grounds it may be expected that the spin rotation 
symmetry U{1) in the xy plane is restored in one space 
dimension, i.e., (si*) — (s 2 %) — 0£ This scenario is sup- 
ported by the renormalization-group analysis of similar 
phases in quantum rotor models^ implying that the true 
tranverse long-range magnetic order in the classical C 



phase is transformed to a quasi-long-range xy order in the 
quantum system. On the other hand, the spin-stiffness 
constant p s remains finite in both magnetic phases as 
well as at the transition point a c \. Following the ter- 
minology of Ref. |5j, the quantum C phase may be called 
partially polarized ferromagnetic phase, as the ferromag- 
netic moment Mq is less than the maximal value s\ — s 2 
in the ferrimagnetic phase (see Fig- El) - This quantum 
state may also be classified as a kind of ferromagnetic 
Luttinger liquid,- We postpone the detail analysis of this 
quantum phase for future studies. 



III. QUANTUM PARAMAGNETIC PHASE 

Now, let us turn to the region a > a c2 of the phase 
diagram characterized by St — 0. It is instructive to 
rewrite Q in the following form 



H = ^ (JiS n • S n+ i + Jj_si„ • s 2n ) 

71=1 

L 

n ' s ln+l + S 2 „ • S 2 „ + i) . (3) 



The operator 2si„-s 2n = (S„) - si(si + 1) - s 2 (s 2 + 1) is 
a conserved quantity for J 2 = 3\ and at this special point 
the low-energy physics of Q is described by the antiferro- 
magnetic spin-S' Heisenberg chain: (S„) 2 = 5(5+1), 5 = 
si — s 2 , Sx — S2 + 1, . . . , si + s 2 . In the case of special in- 
terest (si, s 2 ) = (1, i), and for relatively small interchain 
couplings Ji J$ 1.59 Ji, we have numerically found that 
all the rung spins are characterized by 5 = |, so that in 
the low-energy sector the ladder model (JJJ is equivalent 
to the 5=| antiferromagnetic Heisenberg chain. 

The above statements concern the special point J 2 = 
Ji where 5 is a good quantum number. As an example, in 
Fig. |21 we show the energies of the lowest excited states 
(L = 12 and 14) for a = 0.4 and 1. Apart from the 
k = state, which is characterized by the total spin 2, the 
lowest excited states are triplets above the singlet ground 
state. This structure of the low-energy spectrum is valid 
in the whole region a > a c2 up to the limit a = 00 where 
the system is composed of two independent si = 1 and 
s 2 = i antiferromagnetic Heisenberg chains. In accord 
to the generalized Lieb-Schultz-Mattis (LSM) theorem, 11 
it is natural to suppose that the gapless linear structure 
of the spectrum around k = 7r survives away from the 
point J 2 = J\. 

To study the properties of the quantum paramagnetic 
phase in the whole region a > a c2l we may compare 
the finite-size scaling properties of the ground state and 
the lowest excited states with those based on the SU (2) 
Wess-Zumino-Witten (WZW) nonlinear a model. This 
model with the topological coupling ko = 1 is believed 
to describe the antiferromagnetic Heisenberg chains with 
half-integer site spinsJ^.. In the following we restrict our 
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FIG. 3: ED results for the lowest excited states of the sys- 
tem (si,S2) = (1, \) for two different frustration parameters 
(a = 1 and 0.4): L = 12 (open circles), L = 14 (filled cir- 
cles). Curves represent the two branches of spin- wave modes 
in the paramagnetic Brillouin zone, as obtained from the lin- 
ear spin-wave theory. The spin-wave energies are multiplied 
by the normalization factor v s /v aw — 1.29 where v s = 3.87 
and Vs-w = 2S = 3 are the spin-wave velocities, as obtained 
by the density matrix renormalization group method— and 
the linear spin-wave theory. Note that the classical transition 
point at a — 0.4606 appears as an instability point for the 
lower spin-wave branch. 

analysis to the (s±, S2) = (1, \) system. According to the 
conformal field theory, the ground state energy Eq (L) of 
a periodic system with length L is given by the following 
expressioi: 

T 



10 



6L 2 



1 



0(g 4 ) 



L 4 



(4) 



Here £0 is the ground state energy per rung in the thermo- 
dynamic limit, v s is the spin- wave velocity, and g = g(L) 
is the effective coupling constant of the marginally irrel- 
evant operator — 2ng3L • Jr at the length scale L. Jl 
and Jfl are the conserved current operators for the left 
and right movers in the WZW theory. The L -4 contribu- 
tion comes from irrelevant operators. The coupling g is 
defined by the renormalization-group (RG) equationi&i^ 



1 

9 



1 



In g = In ■ 



L 



(5) 



where L c is a non- universal effective length scale depend- 
ing on the microscopic model. An iterative solution of (0 
yields the following expansion for g 



9 



1 Inln(Z/L c ) 
HL/L C ) ~ 2\n 2 {L/L c ) 



O 



(6) 



so that the marginally irrelevant operator introduces log- 
arithmic corrections in Eq. if4"jl. 

The energy of the lowest triplet excitation E t (L) with 
momentum k = 7r can be expressed in the form 



Et — Eq 

L 



2itv s 



+ b ig 2 + 0(g 3 ) 



L 4 



(7) 



At moderate length scales (L=8,10,12, and 14), the cou- 
pling constants g(L) in (0J and JJJ may have different 
values, so that instead of L c we introduce the effective 
length scales Lq and L t for the ground state energy and 
the energy of triplet excitations. 
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FIG. 4: Reduced energy gap between the lowest excited 
state at n — 2n/L and the ground state. The interpola- 
tion of ED data is performed by the ansatz E(2n / L) /2ttL = 
I]^ =1 c„L~ 2n . For a = 0.45,0.5, and 1 this yields, respec- 
tively, the spin-wave velocities v s = 1.19, 1, 78, and 3.81. Note 
that the interpolation function for a = 0.4 has a positive cur- 
varute. 

The parameter v s in Eqs. (QJ and Q can be inde- 
pendently determined from the scaling of the reduced 
gap E(2tt/L)/2ttL, where E(2tt/L) is the energy gap be- 
tween the lowest excited state with momentum tt — 2it/L 
and the ground state. Using the interpolation ansatz 
E(2ir/ L) /2ttL — J2^i=i c nE~ 2n , we find, in particular, 
the estimate v s = ci = 3.81 at a = 1. This is close to the 
density matrix RG result 3.87 ± 0.02 for the antiferro- 
magnetic S = § Heisenberg chain. 12 In Fig.0]we present 
the interpolation curves for different values of the frus- 
tration parameter a. Excluding the point a — 0.4, our 
estimates for v s can be well interpolated by the ansatz 

v s =Vi(a - a 2c y + v 2 {a - a 2c ) +0[(a — a 2c ) 2 ] (8) 



|. The linear spin- wave 



up to a = 1, provided that 7 w 

theory gives the exponent 7 = i. Note that the spin- 
wave ansatz (JSJ assumes that the velocity v s vanishes at 
the critical point a c2 . Of course, the above interpolation 
of ED data can not definitely confirm such an assump- 
tion, although the apparent change in the curvature of 
E(2tt/L)/2ttL vs 1/L 2 close to a = 0.4 gives some indi- 
cation in favor of v s (a c2 ) = 0^ 

Having the parameter v s for different a, now we can 
interpolate the ED data for E (L) and E t (L) by using the 
scaling expressions (0J and (JJJ. The fitting parameters 
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FIG. 5: Scaling of the ground state energy Eo(L)/L at a = 
0.45 and 0.5. The interpolation of ED data is performed by 
Eqs. and 0. 

in are e o, a i an d Lq. Alternatively, in Eq. (7J the 
fitting parameters are &i, &2, and if. As an example, in 
Fig. [S] we present the interpolation curves Eq(L)/L vs 
1/L 2 for frustration parameters a — 0.45 and 0.5. Using 
the RG equation JSJ for g(L), the best fit at a = 1 is 
obtained for e = -2.3290, Lq = 0.94 and a x = -20.5. 
This is in accord with the density matrix RG result eo = 
— 2.32833. 12 « The parameter Lq = 0.94 corresponds to an 
effective coupling constant 5(10) = 0.35. Performing the 
fits down to a = 0.4, we observe that the characteristic 
length Lo remains almost unchanged (excluding the point 
a = 0.4 where formally Lq — > oo). An interpolation 
procedure using only the leading term in the logarithmic 
expansion |JB} produces similar results, although with a 
slightly larger effective length Lq (= 1.08). 
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states E S (L) may belong to different conformal towers. 
That is why, instead of utilizing the combination of E t (L) 
and E S (L) which eliminates the O(g) correction^ we 
have performed the interpolation directly with Eq. (J7J 
by using the logarithmic expansion © up to second order 
in 1/ ln(L/Lf). The results presented in Fig. Eland Table 

TABLE I: Values of the coefficients in (@J and Q giving the 
best fit to the ED data. The values of v s correspond to the 
ansatz ® with 7 = § giving the best fit to the estimates from 
Fig- B] Lq = 0.94 for all values of a presented in the Table. 



a 


Vs 


eo 


ai 


L t 


bi 


b 2 


0.42 


0.71 


-1.5608 


-26.8 


2.20 


-2.2 


318 


0.45 


1.21 


-1.5873 


-26.0 


2.20 


-1.8 


449 


0.50 


1.77 


-1.6373 


-25.5 


2.12 


-1.5 


495 


1.00 


3.87 


-2.3290 


-20.5 


2.06 


-1.4 


942 



0.001 



0.008 
1/V 



0.012 



0.016 



HI imply that the effective coupling at a given length scale 
g(L) exhibits only a small increase^ when approaching 
the critical point a C 2, in agreement with the interpola- 
tion result for Eq(L). On the other hand, as a — ► a C 2 
one indicates a monotonous growth of the 0(g 2 ) contri- 
bution to the scaling dimension. The scaling behavior of 
Eq(L) and E t (L) for a > 1 qualitatively reveals the same 
properties. 



IV. CONCLUSIONS 

In conclusion, we have examined the impact of mag- 
netic frustration on the ground-state phase diagram of 
two coupled mixed-spin (si,S2) = (U f ) ferrimagnetic 
Heisenberg chains. The analysis of ED data implies an 
interesting phase diagram containing the ferrimagnetic 
phase and a singlet paramagnetic phase exhibiting the 
characteristics of the critical spin-liquid phase in half- 
integer-spin antiferromagnetic Heisenberg chains. Both 
phases are separated by a tiny but finite region occu- 
pied by a critical partially-polarized ferromagnetic phase. 
It is natural to expect similar phase diagrams for the 
whole class of frustrated (si,S2) two-leg ladders with 
half-integer rung spins. 



FIG. 6: Scaling of the triplet gap A t /L = [E t (L) - E (L)]/L 
for the frustration parameters a — 0.42, 0.45, 0.5, and 1. The 
fit of ED data is performed with the scaling formula Q and 
by using the logarithmic expansion © up to second order in 
l/In(L/Lt). 

The O(g) correction to the scaling dimension in @ 
makes the fit of our ED data more intricate. Moreover, 
we have found that for L < 14 the lowest singlet excited 
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